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We study the Schwinger effect, in which the external field having a spatiotemporal profile creates 
electron-positron pairs via multidimensional quantum tunneling. Our treatment is based on the 
trace formula for the QED effective action, whose imaginary part is represented by a sum over 
complex worldline solutions. The worldlines are multiperiodic, and the periods of motion collectively 
depend on the strength of spatial and temporal inhomogeneity. We argue that the classical action 
that leads to the correct tunneling amplitude must take into account both the full period, T and 
the first fundamental period, Ti. In view of this argument we investigate pair production in an 
exponentially damped sinusoidal field and find that the initial momenta for multiperiodic trajectories 
lie on parabolic curves, such that on each curve the ratio T/Ti stays uniform. Evaluation of the 
tunneling amplitude using these trajectories shows that vacuum decay rate is reduced by an order 
of magnitude, with respect to the purely time-dependent case, due to the presence of magnetic field. 


The surge of interest on the Schwinger effect, the non- 
perturbative production of electron-positron pairs from 
vacuum in an external electric field, has yielded new 
insights into this peculiar yet unobserved prediction of 
QED [U [2]. Going beyond uniform field approximation, 
computations of vacuum decay rate show that mean num¬ 
ber of produced particle pairs depends nontrivially on 
the shape, cycle structure and polarization of the exter¬ 
nal pulse [3Hzl- While such investigations predominantly 
deal with the time-dependent electric fields, laboratory 
fields are usually composed of Gaussian or x-ray beams, 
which have spatiotemporal profile and include magnetic 
fields as well. This makes the formulation of vacuum de¬ 
cay in multidimensional electromagnetic fields essential 
for the fully realistic treatment of the problem. 

The technical challenge is to compute the imaginary 
part of QED effective action, Imr[74^], which requires 
the knowledge of vacuum persistence amplitude in the 
background gauge field, A^{x). The standard approach 
is to use the S-matrix formalism, which relates Imr[74^] 
to the WKB coefficients of the vacuum state by a Bogoli- 
ubov transformation [HI [9]. While this approach works 
perfectly well for one dimensional external fields, its ex¬ 
tension to higher dimensional, nonseparable backgrounds 
is exceedingly difficult. One way to circumvent this diffi¬ 
culty is the inverse scattering method, where one starts 
with a plausible ansatz for the Dirac equation to estab¬ 
lish a physical set of gauge configurations pH] . A more 
direct approach involves brute force integration of Dirac 
equation on a spatial grid. An application of this to 
Schwinger effect was given in m- In connection with 
the standard scattering methods, works aiming at gener¬ 
alization of quantum kinetic equation to 1+1 dimensions 
can be found in [12] . Apart from the technical subtleties 
of the chosen method, the integrability of multidimen¬ 
sional Dirac Hamiltonian becomes an important aspect 
of the problem, because it relates to the question of what 
the conserved quantities are in a given background. Ap¬ 
parent translational symmetry of the Hamiltonian in one 
dimensional backgrounds makes the quantization of the 


created particle pairs by the conserved momentum (or en¬ 
ergy) straightforward. In the multidimensional setting, 
identification of particle content remains elusive within 
the framework of Bogoliubov-type transformations. Evi¬ 
dently, these problems also appear in semiclassical WKB 
analysis, which uses analytic continuation of WKB solu¬ 
tions in the complex domain [T3l [T^ . 

Here, we give semiclassical treatment of the multidi¬ 
mensional vacuum pair production by using the world¬ 
line formulation of QED [TSHU]- The worldline approach 
could be considered more advantageous with respect to 
conventional WKB methods for two reasons. First rea¬ 
son is that in the worldline language Imr[A^] is repre¬ 
sented by a path integral over closed trajectories in space- 
time, thus the formalism admits a natural multidimen¬ 
sional description. Secondly, no specific choice of ansatz 
and Bogoliubov transformation are needed; calculation of 
Imr[A^] is relegated to finding periodic, tunneling tra¬ 
jectories, which are also referred as worldline instantons. 
Basic formalism is not just relevant to Schwinger effect, 
but it also has the potential to deal with multidimen¬ 
sional tunneling problems in general. 

Main ingredient of the following analysis is the QED 
analog of widely used trace formula for the Green 
function p!8l [19]. The trace formula can be obtained 
upon performing saddle point approximation to r[A^], 
which can be written as a sum over closed orbits = 

n = c=l): 


( 1 ) 

The matrix A comes from the saddle point expansion and 
corresponds to the density of the trajectories which start 
at the same point with different initial momenta. For 
each periodic trajectory labeled by p, the Morse index, 
rup is given by the number of negative eigenvalues of the 
determinant, whereas the spinor term contributes a +/— 


j ^^ p-iW[Tp]-imp7T/2 

^-tr 


\/det A 




pTp 

W[Tp] = / p^x^^du, 
Jo 


= 2 [ 7 ^ 7 "] 






2 


sign. The determinant can analytically be obtained when 
equations of motion separate, and the trajectory admits a 
simple form. Here, we will work with nonseparable cases 
and focus our attention to the classical action, H^[T], 
which is sometimes referred as the Hamilton’s character¬ 
istic function. Semiclassical value of H^[T] is given by 
the classical tunneling trajectories, which are imaginary 
proper time solutions of the force equation: 

= i\e\F^jy{x)x'", u^iu ( 2 ) 

Periodic solutions of © also exist for real u [20]. Com¬ 
posite worldlines which are made of x{u) and x{iu) are 
responsible for possible interference effects, but their con¬ 
tribution to pair production rate is controlled by the tun¬ 
neling segments x{iu). For this reason, we are interested 
only on imaginary proper time solutions, as they are ex¬ 
pected to yield a pretty accurate estimate in the non- 
perturbative domain. In the following, we argue that 
worldline trajectories are multiperiodic in spatiotemporal 
backgrounds. Starting with a simple model, we discuss 
the implications of bounded motion in the multidimen¬ 
sional setting, and move onto a more realistic scenario 
incorporating the magnetic field. 
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FIG. 1. Worldline trajectories in the background field it , as 
a function of propertime u. Temporal parameters are fixed 
as Fo = 0.1m, uj = 0.1m. Spatial parameters are chosen 
as Fq = 0 (top-left), Fq = 0.01 m^, k — 0.001m (top-right) 
Fq = 0.01 m^, k = 0.01m (bottom-left), and Eq = 0.03m^, 
k = 0.01m (bottom-right) 



We first consider the combination of two pulse-shaped 
electric fields. Total field along x^ is given as: 

F E' 

^xoixs) = —tanhwxo, A^^sixo) = -^tanhkxs 

UJ K 

Fa ,3 (xo, xs) = Eq sech^ ujxq + Fq sech^ kxs (3) 


where k and uj represent spatial and temporal width re¬ 
spectively. To illustrate the effect of spatial inhomogene¬ 
ity on the closed orbits, we keep the temporal parame¬ 
ters fixed, and vary Fq and k. As observed in [20], in the 
purely time-dependent case (Fq = 0) tunneling trajecto¬ 
ries are periodic with the classical period being T, and 


quantized by the canonical momentum ps. Introduction 
of spatial inhomogeneity causes periodic trajectories to 
have a second oscillation period T, which envelopes the 
oscillations with a smaller period, Ti (Fig. . These 
trajectories are quasiperiodic and form invariant tori in 
phase space. Given the time component of the above field 
is dominant (Fq < Fq), such bounded trajectories gen¬ 
erally persist for a finite range of momenta, when k < uj. 
Our key observation is that in the constant spatial field 
limit: /c ^ 0, and also in the limit: Fq ^ 0, T in¬ 
creases and effectively goes to infinity, leaving a single 
period Ti = T. This indicates finite values of Ti and 
T genuinely depend on the interplay between the tem¬ 
poral and spatial adiabaticity parameters of the external 
field. This applies not just to this particular case but 
holds general validity. Few other examples include pulse 
configurations such as: 


Ex 


Ex 


= FQe-^'"’3-^'^o^ =FQe-^'"’3sechcJXo^ 
= Fq -h F’o 


Note that Hamiltonian, EL = ll2 — ejc A^{x))^^ for 
the external field in has no longer translational 
symmetry along x^. On the other hand, existence of 
quasiperiodic motion tells us that system possesses dy¬ 
namical (hidden) symmetry. This is because Poincare 
sections of the quasiperiodic orbits form closed curves in 
the reduced phase space, and imply the existence of an 
integral of motion C, in addition to EL. One may hope 
in this case the corresponding quantum system possesses 
simultaneous eigenstates of constants of motion by the 
virtue of Liouville integrability. Quantum tunneling via 
multiperiodic motion can then be interpreted as the cre¬ 
ation of particle states quantized by EL and C. 

The tunneling amplitude for the vacuum decay is gov¬ 
erned by the action, whose evaluation on multiperiodic 
solutions must be handled with care. Direct approach 
would involve computing IF[T], where T simply gets re¬ 
placed by T. This however leads to an undesirable lim¬ 
iting behavior for the tunneling amplitude. To explain 
we reconsider the toy model above. For a weak spa¬ 
tial inhomogeneity (Fq ^ Fq) we have the second pe¬ 
riod extended over a very large proper time scale com¬ 
pared to Ti, so we have Ti ^ T. In this regime, the 
smaller period almost coincides with the period of purely 
time-dependent background, whose action is denoted by 
W[T]. On the other hand, the evaluation of W[f] yields 
^ T/TiW[T]. Thus the corresponding vacuum decay 
rate is many orders of magnitude smaller than 
which gives the decay rate for the time-dependent field. 
Such a huge difference in tunneling probabilities is quite 
unexpected for a weak spatial perturbation. In fact, in 
the limit: Fq ^ 0, T gets infinitely large and so W[f] will 
get infinitely large. This basically leads to vanishing pair 
production probability, which is clearly in contradiction 
with the purely time-dependent case. To remedy this 
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FIG. 2. External fields represented by 0. Field parameters 
are given in the text. 


situation tunneling treatment of the trajectories should 
take the internal oscillation period Ti into account. To 
incorporate Ti into this picture, we rewrite the effective 
action as: 


r[^] = 2*^ Log 


+ Ip 

dTi 

/ f r exp 

Jo T, 


-^(m2 + ^')T 
2T 




]p=r{d^ - ieA^), W[T^, T] = F p^x^^du (4) 

1 Jo 


Integration variable in the definition of the logarithm 
above was chosen to be Ti, whose appearance would pro¬ 
duce the desired limiting behavior for the tunneling am¬ 
plitude. Through redefinition of the exponent in second 
line ofand fixing the upper limit of the ordering param¬ 
eter u as T, the trace operation is to be performed over 
the space of multiperiodic trajectories, for which Ti/T 
is a constant factor. With this in mind and making the 
change of variable (Ti T) in the integration, saddle 
point approximation over T can be carried out in the 
usual way. At first look the transformation of the expo¬ 
nent inj^is very similar to gauge fixing of the effective ac¬ 
tion; but it does not correspond to a simple deformation 
of the period, T, rather it represents IT[Ti, T] as an av¬ 
erage over the internal cycles of the trajectory. Note that 
the evaluation of kF[Ti, T] on quasiperiodic solutions not 
so surprisingly yields a small, unphysical imaginary part, 
because the trajectory does not close back on itself at 
some finite T. Yet this imaginary part tends to vanish in 
the limit T ^ oo, which precisely corresponds to taking 
the average of a quasiperiodic function. The question of 
whether such quasiperiodic averaging may contribute, or 
assist to quantum tunneling might be interesting in its 
own right. In the remainder of this work however, we 
will focus on the closed orbits. 

In the following we introduce spatial variation to the 
electric field: Ei[xq) = coscj^o, whose corre¬ 

sponding vacuum decay rate was analyzed in [5]. This 


field configuration approximates the experimental setup 
where two counter-propagating short laser beams form 
a standing wave. Interaction region is conveniently cho¬ 
sen along the beam, on the plane where magnetic fields 
cancel. Here, we release this restriction by taking into 
account the spatial profile of the external field, along the 
beam direction, x^. For this, we consider the transverse 
field: 


X 3 ) 

Exii^o, X 3 ) 
Bx^ixo, X 3 ) 

f(xo) 


-Eoe 2/2 ^ 2 (e f{xo) + c.c) , 

2 2 
_ ^0 _ ^3 

Eqc ^ ^ COS (uj xo — k xs), 
Eoe~^~^ g{x 3 )f{xo) + c.c) , 


0FTErf 


gixs) = ^+ ik (5) 


2^2 


with the frequency, w and the wavenumber, k (Fig. |^. 
The temporal and spatial width are respectively given by 
r and a. The exponential terms above make the external 
field finite, whereas oscillatory terms simulate the inten¬ 
sity variations typically seen in the Gaussian or x-ray 
pulses. Note that above field introduces a source current 
along xi. This current may in principle affect the tun¬ 
neling amplitude, because it contributes to the evolution 
of the momentum operator, dxi/du. However, the effect 
of the source term remains negligible, and classical so¬ 
lutions are expected to dominate tunneling, as long as 
the external field extends over a distance larger than the 
Compton wavelength [ 21 ]. 

The existence of invariant tori now depends on the 
relative magnitude of both a and k with respect to the 
temporal parameters. Quasiperiodic trajectories gener¬ 
ally occur in the parameter region where ka < ujt and 
k < uj. Here, we look for the closed orbits when the elec¬ 
tric field is maximum, so we fix the initial positions as 
xo( 0 ) = ^ 3 ( 0 ) = 0 , and vary the conserved momentum 
xi[0] = —ipi, and the initial velocity X3(0) = —ipo- The 
remaining initial condition for ± 0 ( 0 ) = —ipo is fixed by 
the constraint: Xq{u) —x\{u) —x\{u) = —m^, where m is 
the electron mass. For the sake of comparison, we specify 
the values of ^ and cj, in accordance with [5] . In terms 
of the normalized mass (m = 1 ) we fix the field param¬ 
eters as: Eq = 0.1 m^,T = 100 uj = 0.03m. To see 
the full effect of spatial dependence and the arising mag¬ 
netic field, we will work in the parameter region where 
temporal and spatial inhomogeneities become equally im¬ 
portant. We accordingly choose the spatial frequency to 
he k = 0.01m and a = 100 m“^. Multiperiodic orbits are 
located by making use of a search algorithm that scans 
through the region: pi G ( 0 ,m),p 3 G ( 0 ,m), until the 
trajectory closes on itself within an accuracy of 10 “^. 
The greater accuracy makes the location of the orbits 
more precise but, from a practical point of view, increas¬ 
ing the accuracy further does not have appreciable effect 
on the tunneling rate. Figureshows the locations of the 
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closed orbits on the momentum plane; multiperiodic tra¬ 
jectories are not isolated, but form a family. Along each 
parabolic curve representing the orbit family, trajecto¬ 
ries grow in amplitude and get steeper until the invariant 
tori break and the motion becomes unbounded. For ev¬ 
ery closed trajectory on the first parabolic curve we have 
T/Ti = 12, on the second, period doubling occurs and 
ratio becomes 23, and on the third curve it becomes 11. 
On the first and third curve T matches with the second 
fundamental period, T 2 , whereas on the second curve we 
have T 2 = T/2. Accordingly, the fundamental frequen¬ 
cies of the system consists of set of co-prime integers, 
which are respectively given as: (12,1), (23,2), (11,1). 
These sets of integers characterize the orbit topology [22]. 
Consequently, 12 is considered as a topological sum, where 
each curve family is labeled by its topological index. 

Before we finally evaluate lF[Ti, T], we analyze the 
spinor term, which leads to an interesting geometric in¬ 
terpretation for the pair production process. Straightfor¬ 
ward calculation of the spinor trace for the field configu¬ 
ration (§ gives (u iu) 


tr 




COS 

T 


2T 


52 


£ = \e\ [Exi du, B=\e\ f du 
Jo Jo 


(6) 


The exponent in has the form of a Lorentz transforma¬ 
tion: as the spin precesses under the influence of external 
field, the amplitude goes under a rotation in 

the complex plane given by the angle Ti/4T j^F^ydu. 



FIG. 3. Location of the quasiperiodic trajectories on pi — 
P 3 plane (shaded). Multiperiodic trajectories are located on 
the curves (top). The figure below shows a typical closed 
trajectory with pi = 0.35351 and pa = 0.42802 


It can readily be shown by integration that for closed 
trajectories the value of Vf ^ is precisely 27rT/Ti. 

Geometric meaning of this result becomes more trans¬ 
parent if one considers the relation between \/f ^ 
and the curve invariants of the motion with the aid of 
Frenet-Serret formulas. For instance, in the case of sin¬ 
gle dimensional inhomogeneities geodesic curvature of a 
closed orbit is given by \e\Ex^- Consequently, the ar¬ 
gument £ is nothing but the total curvature, /C, whose 
value is 27r, regardless of the shape of the electric held. 
This is because instanton trajectories with 1 degree of 
freedom are simple closed loops, which can be continu¬ 
ously deformed to a circle. With 2 degrees of freedom 
trajectories become non-planar, and geodesic curvature 
can not be solely given in terms of the held strength. 
Also, total torsion along the closed orbit may not vanish, 
unless the given trajectory is symmetric under reflections 
with respect to the normal plane at = 0. In the general 
case the relation between the held strength and the curve 
invariants can be given by: 

\\e^F^^du\\^^/K^-T^ (7) 

Jo 

The norm above gives the total rotation angle in inho¬ 
mogeneous fields, whereas the right-hand side could be 
regarded as the total invariant curvature of the trajec¬ 
tory. For the closed curves considered here 1C is very 
close to 2'kT jTx^ whereas the total torsion, T, vanishes. 
For slowly varying transverse electromagnetic fields, it 
could be said without loss of generality that the imagi¬ 
nary part of T[A^] follows from the quantization of the 
rotation angle by the integer multiples of 27r. But unlike 
the uniform or purely time-dependent cases, quantization 
condition in spatiotemporal backgrounds depends on the 
set of values that T/Ti can take. If the semiclassical 
approximation holds well, it is then evident that the ex¬ 
istence of worldline trajectories with quantized invariant 
curvature is intimately connected to the nonperturbative 
pair production. 

Returning back to evaluation of LF[Ti, T], we see that 
for a single closure the spinor trace yields cos tt, which can 
be absorbed into the prefactor. This further simplifies 
imaginary part of the effective action to: 

Inir[A^] « ^ W[T,, f] > 0 (8) 

P 

Having obtained the closed orbits, decay rate can now 
be determined upon evaluation of the action. For con¬ 
venience, we use the configuration space where LF[Ti, T] 
reads 

W[Ti, f] = -Ti/f [ {xl - ±1) du. (9) 

Jo 

Figure [^ shows vacuum decay rate, i.e tunneling proba¬ 
bility for each family of closed trajectories. Decay rate is 
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FIG. 4. Tunneling rates for each family of closed trajectories 


consistently higher for the family of orbits with relatively 
lower momenta, but it surprisingly increases along all the 
curves, until the invariant tori break. This behavior sug¬ 
gests that decay rate is maximized on the boundary of 
tori and it would be interesting to see whether this re¬ 
mains true for other cases. Comparing the results with 
purely time-dependent counterpart of ([^, we see that 
pair production rate for the overlapping range of canon¬ 
ical momentum drops down by an order of magnitude 
[5]. The fall of vacuum decay rate in the presence of a 
perpendicular magnetic field can in fact be seen directly 
from the pair production formula for the uniform field. 
By using the peak values of electric and magnetic field in 
the invariant field strength, it is easy to see that uniform 
field approximation also yields around an order of mag¬ 
nitude drop for the decay rate, with respect to the purely 
electric field case. This result can intuitively be under¬ 
stood by using the qualitative description of the virtual 
particles: the external electric field along xi separates 
the virtual electron-positron pairs, whereas the magnetic 
field applied in the perpendicular direction tries to bring 
these virtual pairs together. 

To summarize, worldline method relies on the exis¬ 
tence of bounded motion. In this respect, integrability 
of the Hamiltonian plays an important role in the pair 
production process. Given that H is integrable, the use 
of multiperiodic trajectories in quantum tunneling natu¬ 
rally leads to the quantization of the fermionic modes by 
ii and C. For the external field given in ([^, multiperi¬ 
odic trajectories form a one parameter family due to the 
existence of an additional constant of motion, pi. But it 
is useful to keep in mind that tunneling may also occur 
via isolated trajectories. To make a final remark on the 
use of transformation in Q, we would like to point out 
that averaging methods similar in spirit have been used 
in atomic molecular physics to obtain the spectrum of 
multidimensional bound systems. For instance, the eval¬ 
uation of action by using the caustics or the Poincare sec¬ 
tions of quasiperiodic trajectories has successfully gener¬ 
ated the energy eigenvalues of multidimensional systems, 
such as 2D coupled harmonic oscillators [23] . Here in our 
approach averaging arises naturally, with only requiring 
the periods of motion. 

The worldline picture provides a valuable semiclassical 
tool to calculate the nonperturbative decay amplitudes 


yet there are several noteworthy aspects of the problem 
that remain unaddressed. The first one is the prefac¬ 
tor contribution to the decay amplitude. A numerical 
method to obtain the prefactor and Maslov index in one 
dimensional inhomogeneities was developed in [16]. Ex¬ 
tension of such method to spatiotemporal backgrounds is 
necessary for the completeness of the worldline approach 
and will be the subject of a future work. Another issue 
is related to the constant of motion, C. The determina¬ 
tion of C from the given orbit data would be a desirable 
asset for the worldline method because the values that 
C takes, in conjunction with the topological index, can 
reveal a good deal of information about the spectrum of 
the created particles. Final and perhaps more interest¬ 
ing aspect is the quantum interference effects seen in the 
particle spectrum. Basic mechanism behind the interfer¬ 
ence phenomenon, in time-dependent fields for instance, 
can be understood in terms of worldlines that are glued 
together at critical (conjugate) points on the complex 
time plane. Multiperiodic counterparts of such compos¬ 
ite orbits are expected to come with a rich topological 
structure, but whether this basic procedure can be ex¬ 
tended to the multidimensional cases remains to be seen. 
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